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We extend the pseudoparticle transport description of the Hubbard chain to all energy scales. 
In particular we compute the mean value of the electric current transported by any Bethe-ansatz 
state and the transport masses of the charge carriers. We present numerical results for the optical 
conductivity of the model at half-filling for values of U/t = 3 and 4. We show that these are in 
good agreement with the pseudoparticle description of the finite-energy transitions involving new 
pseudoparticle energy bands. 
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I. INTRODUCTION 

The electronic transport properties of strongly cor- 
related electron systems in low-dimensional conductors 
have been of experimental and theoretical interest for a 
long time. Low dimensional conductors show large de- 
viations from the usual-single particle description which 
suggests that electronic correlations might play an im- 
portant role Q. Although the Hubbard chain has been 
diagonalized by Lieb and Wu long ago, the involved 
form of the Bethe ansatz (BA) wave function has pre- 
vented the calculation of dynamic response functions, 
like the optical conductivity a(u>). Several approaches, 
using perturbation theory || , the pseudoparticle formal- 
ism m\, bosonization f|, scaling methods 0, and spin 
wave theory [Q have been used to investigate the low en- 
ergy behavior of <j(uj) away from half-filling and at the 
metal-insulator transition ^ . Interesting information on 
of a{uj) at finite values u> has been obtained numerically 

Recently, it has been possible to obtain analytical re- 
sults for the finite-a; behavior of cr(u>) using a pseudopar- 
ticle description for all the Hamiltonian eigenstates of 
the model [jyj. This theory generalizes, for all energy 
scales, previous results by Carmelo et al. JO], by intro- 
ducing new branches of pseudoparticles. The latter are 
called heavy-pseudoparticles because the eigenstates de- 
scribed by them have an energy gap relatively to the 
ground state. As in the case of the low-energy properties 
of the ID Hubbard model [O], it is possible to write the 
Hamiltonian in terms of a set of anti-commuting opera- 
tors, called pseudoparticle operators, and to generate all 
the model eigenstates from the 5 l O(4) ground state [ fijl . 
We show below that the finite-w properties of the optical 
conductivity at half filling are determined by transitions 
involving the bands populated by the usual pseudoparti- 
cles (called c and s pseudoparticles [|l2|]) and the heavy- 
pseudoparticle bands (called c, 7 and s, 7 pseudoparticles 
0)- 



The paper is organized as follows. In Section II we 
study, within the pseudoparticle formalism, the response 
of the Hubbard ring to a small electric field and com- 
pute the mean value of the electronic current for any 
BA state and the transport masses of the charge carri- 
ers. In Section III we present numerical results for the 
optical conductivity of the half-filled ID Hubbard model 
and compare them with the pseudoparticle-theory pre- 
dictions. In Section IV we summarize our results and 
discuss possible future developments. 



II. RESPONSE TO AN ELECTRIC FIELD: 
PSEUDOPARTICLE DESCRIPTION 



The linear response of the ID Hubbard model to a 
small spatially uniform, time-dependent, electric field 
parallel to the chain is achieved by means of the usual 
Peierls phase substitution of the hopping integral (t — > 
l e ±%<s>(t)/N a y This means that one threads the ring by a 
flux $(i) [HQ. The obtained Hubbard-chain Hamilto- 
nian (with N electrons) in a magnetic field H and with 
chemical potential /1 is written as 

H = -tY / [e +im/N "clc J+la + h, 

3><r 

+ U^D 3 - /i(A a - ^ ™j>) - ™ 3 > , (1) 
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where Dj = [hj^ — 1/2] [r^j, — 1/2], Cj a (cj a ) creates (an- 
nihilates) an electron with spin a (a =|, J, when used as 



operator index or a = ±1 otherwise), h 



3,<? 



- J 



^ja^3" 



is the 



number operator at site j, N a is the number of sites of 
the ring, and CN a +ia = ci CT (we are using units such that 
% = 1, the lattice spacing a = 1, and the electron charge 
e = —1). The form of the interaction term accounts for 
the particle-hole symmetry of the model at half filling 

Although the pseudoparticle description refers to all 
Hamiltonian eigenstates [[TlJ, in this paper we restrict 
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our study to the Hilbert subspace involved in the zero- 
temperature conductivity spectrum. Since the current 
operator commutes with the eta-spin and spin operators 
r) 2 , fj z and S 2 , S z , respectively, the Hilbert subspace is 
in the present parameter space spanned by the lowest- 
weight states (LWS's) of the eta-spin and spin algebras 

It has been possible to solve the Hamiltonian (jl|) 
by means of the coordinate BA both with twisted and 
toroidal boundary conditions (both approaches give es- 
sentially the same results) |l4||l7]]. One obtains the en- 
ergy spectrum of the model parameterized by a set a 
numbers {kj, S-yj, Rj,j} which are solution of the BA in- 
teraction equations. Following Takahashi's $(t) = for- 
malism H|, the interaction equations with the flux Q(t) 
read 



and N. 



All the LWS's of the model Cfll are deter- 



k,N a = 27r/; + *(t) -EE 



7 =oy=i 



sin(kj)/u — S^ji 



(7 + 1) 



7=1 j'=l 



sin(kj)/u — Rj,j' 



(2) 



2N a ResirT 1 [{R 1 j + ry)u] = 2nl^ n + 2 7 $(i) 



N r . 



j'=i 



sin(kji)/u — Ryj 



7' = 1 j'=l 



(3) 



and 



(7 + 1) 



®7+l,7'+l(&M ~ ^7',i') • ( 4 ) 

7'=0 i'=l 

In Eqs. (g), (|), and (g) we have that 6> = 2tan- 1 (x), 
u = U/At, and the functions 6 7l7 ', and Q 7 +i .y+i have 
the same definition as in Takahashi's paper |lq| . The fol- 
lowing definitions, A" +1 /u = S 7 j (with n + 1 = 7 and 
a = j), AJ l /u — Ryj (with n = 7 and a = j) with 7 

1.2 x for the iV c . 7 sums and 7 = 0, 1, 2, . . . , 00 for 

the N sn sums, allows us to recover Takahashi's formulae 
for $(i) = [p| (often we also use the notation c = c, 0, 
with the c, 7 sums running over 0, 1, 2, . . . , 00). The quan- 
tum numbers J?, ^' 7 , and /j 1 ' 7 are consecutive integers or 
half-odd integers depending on the parity of the numbers 
E 7 =o^7+E 7 =i N c , y , N a -N + N cn , and N-N s , y , re- 
spectively. The conservation of the electrons numbers im- 
poses the sum rules N± = 5Z 7 =i l N c,j + J2 1 =oh + 1 ) N s,j 
and N C = N -2^=! jNc ,7 on the numbers N c , N c ^, 



mined by different occupancies of these quantum num- 
bers. For example, the ground state is described by a 
compact symmetric occupancy around the origin of the 
numbers and Ij' , and by zero occupancy for the num- 
bers i^' 7 and i"J' 7>0 ||. Equivalently, the eigenstates of 
the model can be described in terms of pseudomomen- 
tum {g"' 7 = 27r/"' 7 /A a } distributions, where a = c, s 



and I C A '° = I 



3 j 

The energy eigenvalues are given by 



N c 



N 



e = -2*E c ° s (%) + EE URe v 1 ~ ^(^tj + 

3=1 7=1 3=1 

+ N a (U/4 - M ) + VV(M - (7/2) - /j.oH(N^ - JV ; ) . (5) 

In the limit of a large system (iV a — > 00, N/N a fixed) we 
can develop a generalization for the low-energy Landau- 
liquid description of the ID Hubbard model ]l9|-pT[| and 
of its operational representation p2| for all energy scales 
Within this new general framework we have the 
usual c, = c and s,0 = s pseudoparticles (and their 
energy bands) , whose different occupancies of the pseu- 
domomenta {2irl^° /N a } and {27!"/?' /N a }, respectively, 
describe the gapless low-energy physics. In addition, we 
have the branches (and energy bands) of heavy c, 7 and 
s,7 pseudoparticles, whose different occupancies of the 



pseudomomenta {2ttIj^ >v /N a } and {2tt/ j s ' 7>u /N a }, re- 
spectively, describe excited eigenstates with an energy 
gap relatively to the ground state . 

Combining Eqs. (gj^) there are several interesting 
transport quantities that can be computed. An exam- 
ple of physical interest for the study of metal-insulator 
transitions is the charge stiffness D [gj|,|Q||,[|| . In- 
formation on D is at zero temperature contained in the 
ground state energy Eq through the relation (24] 



D = - 



1 d 2 (E /N a ) 



2d^(t)/N a ) 2 



(6) 



*(t)=o 



This quantity can be computed within the pseudoparticle 
formalism by combining a Boltzmann transport descrip- 
tion for the pseudoparticles with linear response theory 
jl| or by direct evaluation of Eq. (||) H^5|. Another 
important transport quantity is the electric current of 
a given eigenstate, with the operator j reading, for the 
Hubbard model, j — itJ2f=i( c ja c j+ic~ c ]+io c j<?)- Since 
in this model H and j do not commute the BA wave 
function does not diagonalizes simultaneously the Hamil- 
tonian (g) and the current operator j. Therefore, the BA 
alone can only provide information on the mean value of 
the electric current. This is an important quantity for it 
allows us to compute the transport masses of the charge 
carriers of the system. For the Hubbard chain, the carri- 
ers are the pseudoparticles c and c, 7, as it follows from 
Eqs. (&H) and from a alternative Boltzmann transport 



2 



analysis |25j . The coupling constants to the electric field 
are given by C p = 1 and C£ 7 = 2 7 (and C£ 7 = 0). 

The mean value of the current operator < m\j\m > 
transported by the eigenstate \m > of energy E m is given 
by -d(E m /Na)/d(^(t)/N a )\^ t ) =0 . As stated above, the 
LWS's are obtained considering all the possible occupa- 
tion distributions of the pseudomomenta 2ttI°~' 1 /N a . For 
this reason, it is simpler to describe < m\ j\m > in terms 
of the pseudoparticle occupation distributions N a ^(q) 
(q being the pseudomomcntum) . The computation of 
< m\j\m > involves the expansion of Eqs. (§-|^) up to 
first order in the flux <&(i). After lengthy calculations, 
we obtain (in the limit of N a — * oo) 

<m\j\m>= ^2^2 d Q N a, 1 (q)j a , 1 (q) , (7) 

a=c,s 7=0 J—1o,,~t 

where j a<1 (q) = E Q ' =C , S Ey=o C a', 7 ' [v a , 1 {q)5 a ., a >5 1 , 1 > + 
l^J2j=±ijfa,T,a',j'(q,jqF a ' n ')}, and v an (q) are the 
pseudoparticle velocities (defined as de° ai {q) / dq, where 
£ a,-y(q) i s the pseudoparticle dispersion), ja,7;a',7'(<7, q') 
are the pseudoparticle interactions (or /-functions), and 
the momenta qFa,j are the pseudo- Fermi points fy[]. The 
first term of ja.~/(q) is what we would expect for a non 
interacting gas of pseudoparticles and the second term 
describes the dragging effect on a single pseudoparticle 
due to its interactions with all the others. 

We stress that Eq. (0) is valid for all energies because, 
in contrast to the Fermi liquid theory, there is only for- 
ward scattering among the pseudoparticles at all energy 
scales. The charge transport masses m p aa are defined as 
to q,7 = QFarf / [Ca^jotj-yiq)] and contain important physi- 
cal information. An interesting and well known property 
is the metal-insulator transition of the model which oc- 
curs at half- filling and zero temperature It has been 
shown by means of the pseudoparticle formalism [|| that, 
at half-filling, the transport mass of the c pseudoparti- 
cles is infinite (whereas m p sl = 00 always), as it should 
be for an insulator, and that the charge stiffness is also 
zero (D oc 1 /m p Q ) , satisfying Kohn criterion J24|] . 

In contrast to the zero temperature limit, the 0,7 
heavy pseudoparticles play an important role in the 
transport properties of the model. We computed the 
transport masses for the c, 7 heavy pseudoparticles in the 
limit of half filling (in a zero magnetic field) and obtained 
m p 7 — > 00 as n — > 1 . Whether this result can shed some 
light on the relation between integrability and anoma- 
lous transport properties |^6| is something that requires 
future studies. 

III. OPTICAL CONDUCTIVITY AT 
HALF-FILLING 



of the general pseudoparticle theory. At half filling and 
zero temperature the charge stiffness is zero and the 
real part of the optical conductivity Rea{ui) has a fi- 
nite energy gap for all values of U/t. At half filling the 
c, pseudoparticle band is full and excitations can only 
occur between the low-energy a bands and the heavy- 
pscudoparticle bands. Thus, the form of the optical 
conductivity is determined by these type of excitations, 
which are rather unusual. For example, to populate the 
lowest c, 7 energy-band (the c, 1 band) with one c, 1 heavy 
pseudoparticle one has to annihilate two c and one s 
pseudoparticles. The energy involved in this process is 
e°i(0) ~ £0,0(^0,0) - £°c,o(-qFc,o) ~ e% {qFs,o) = 2/x(n = 
1), where [i(n — 1) stands for the chemical potential in 
the limit of half filling which equals the Mott-Hubbard 
gap |pH . In addition, this excitation process is character- 
ized by a topological excitation of the s pseudoparticles 
of momentum —qFs.o @|) ensuring that the overall mo- 
mentum vanishes, as required for the excitations induced 
by the current operator. In general, the process of cre- 
ating one c, 7 heavy pseudoparticle has a minimum en- 
ergy cost of 27/x(l) and a maximum cost of 27^(1) + 8 7 i. 
These transitions are expected to show up in the optical 
conductivity as a set a optical bands. 

In order to investigate this problem, we have studied 
the model numerically JlFJ), by using a quantum Monte 
Carlo technique for systems with 32 sites |27l, which is 
considerably larger than previous exact diagonalization 
studies |^||- The current-current correlation functions 
are computed in imaginary time and continued to real 
frequency using the so called maximum-entropy method 
p8j. The resolution of this method is limited but, never- 
theless, gives a semi-quantitative confirmation of a multi- 
band optical conductivity. 

In Fig. we show our results for U/t — 3 and 4 where 
two peaks are clearly resolved. As indicated in the figure, 
the lower edges of the spectra agree very well with the 
predicted gaps, and the second peak appears above the 
predicted lower edge of the second band, as obtained from 
the pseudoparticle theory, which also allows the evalua- 
tion of the critical exponents associated with the conduc- 
tivity edges JllJ . Note that the weight of the first band is 
concentrated towards the lower end of the allowed band 
of width 8t. The second peak seen in our results is prob- 
ably dominated by the second optical band, but likely 
contains contributions also from the tail of the first band 
(higher bands cannot be resolved due to the limitations 
of the method). 

As discussed in detail in Ref. poj] , for large values of 
the Hubbard interaction U only a single broad peak of 
width « 8t is observed, likely due to the low intensity of 
the higher bands. 



As discussed in the previous Section, some transport 
properties of the model can be understood from the use 
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IV. SUMMARY AND OPEN QUESTIONS 

In this paper we developed a general theory, valid for 
all energy scales, for the electronic transport properties 
of the Hubbard chain (details of the calculations will be 
presented elsewhere j25|). The fact that [H,j] is non 
zero does not allow us to obtain information on the cur- 
rent operator j it self from the BA equations but only on 
its mean value. (The conductivity exponents can be ob- 
tained by combining BA and properties following from 
the model infinite conservation laws.) Still, this is an 
important quantity, since it is related to the transport 
mass of the charge carriers (the c, 7 pseudoparticles, with 
7 = 0, 1, ... , 00) which is believed to be important in the 
transport properties of the model at finite temperatures. 
This question is now an issue of ongoing research. 

In addition, using maximum-entropy analytic continu- 
ation of quantum Monte Carlo data, we computed the op- 
tical conductivity of 32-site Hubbard chains for U/t = 3 
and 4, and found a two-band structure with band edges at 
energies consistent with our pseudoparticle theory. This 
agrees with the conductivity absorption edges obtained 
directly from the pseudoparticle theory. 
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FIG. 1. Optical conductivity results for 32-site Hubbard 
chains at half-filling, for values of U/t = 3 (top) and 4 (bot- 
tom), computed using maximum-entropy analytic continua- 
tion of quantum Monte Carlo data. The solid and dashed 
vertical lines indicate the theoretical predictions for the lower 
edges of the first two bands. 
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